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Исследуется возможность определения полей подвижности воздуха, температур и концентраций вредных 
веществ на основе конечно-элементного моделирования процессов массопереноса заг| рязнений в производ- 
ственной среде с учётом завихрений воздушных ПОТОКОВ. 
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Введение. В работе [1] определение полей концентраций вредных веществ, температуры и по- 
движности воздуха в производственной среде исследуемых помещений осуществлялось на основе 
моделирования тепломассопереноса для потенциального безвихревого поля. Рассчитанные пара- 
метры производственной среды достаточно хорошо согласовывались с экспериментальными дан- 
ными, полученными в помещениях с невысокой плотностью размещения стационарных источни- 
ков загрязнения. 

Однако у источников загрязнения, где, как правило, размещаются рабочие места операто- 
ров и местные вентиляционные отсосы, создаются условия для возникновения турбулентных га- 
зовоздушных потоков, которые существенно изменяют картину распределения подвижности воз- 
духа и загрязнений в помещении, что необходимо учитывать в разрабатываемой модели. В рам- 
ках подхода, изложенного в [1], задача корректно не решается. 

Таким образом, для замкнутых производственных помещений с высокой плотностью раз- 
мещения источников загрязнений в условиях работы активной вентиляции возникает проблема 
точного определения полей подвижности воздуха с учётом вихревых движений, точность опреде- 
ления которых влияет на результаты расчёта полей концентраций и температур. Поэтому необхо- 
дим поиск математических моделей, которые могут описывать исследуемые процессы с высокой 
ТОЧНОСТЬЮ. 

Обзор и анализ отечественных и зарубежных литературных источников [2, 3, 4], посвящённых 
процессам массопереноса веществ в замкнутых средах, показал, что реализация метода «вектор 
завихрённости — функция тока» имеет ряд преимуществ. Так, уравнения предлагаемого метода 
подобны по типу (по математическим свойствам), и их численное решение проще, чем решение 
уравнения Навье — Стокса. 
В работе [4] рассмотрен численный анализ конвекции в прямоугольной области с источ- 
ником тепла в условиях внутреннего массопереноса и внешнего вынужденного течения (1). 
Постановка задачи (1), как и многих других задач вязкой несжимаемой жидкости в пере- 


менных (4,0), обладает следующей особенностью. Граничные условия на твёрдой стенке зада- 


ЮТСЯ ТОЛЬКО ДЛЯ функции тока, а не для вихря, который определён ЛИШЬ внутри области согласно 
0?ф 0? 

2 + 2 
0х? ду 
частности, применяют приближённые граничные условия для вихря. 








уравнению =-@. Для преодоления этой трудности используют различные подходы, в 


10 


Вестник ДТУ. 2012. № 6 (67) 








В данной работе для вихря на твёрдой стенке ставилось условие Тома, которое получа- 
лось из условия прилипания [5]. 


9(и0) 0(ю)_ [5 р ия 








+ 
9х ду 0х? ду? 
0°ф  д?ф 
о =- 
ОХ 0* 
9(иС) ом 2 2 
( и. ОС: +0 (1) 
ох У 55 [0Х? ду? 
И =дф/ ох 
-И=дф/ду, 


гдех,У — координаты Декартовой системы координат; Ш,ИУ — составляющие скорости в проек- 
ции на оси Х,У соответственно; Ф — функция тока; ©@ — завихрённость скорости; С — кон- 
центрация примеси в области решения; О — источник загрязнения. 

Система (1) не может быть решена в общем виде аналитически, для нахождения решения 
необходимо использовать численные методы, реализация которых возможна в программной среде 
НехРОЕ-6.20. 

Постановка задачи. Объект исследования (рис. 1) представляет собой воздушную камеру с 
входными и выходными воротами и расположенным по центру источником выброса оксида угле- 


рода. Скорость газа поступающего в камеру является постоянной в течение всего процесса. 
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Рис. 1. Воздушная камера с источником загрязнения 


Вихревая модель массопереноса вредных веществ. Для универсальности система уравне- 
ний (1) была приведена к безразмерному виду, а безразмерные переменные величины приняли 
следующий вид: 

Х=Х,, У=УШ,, С=С/С,/, 9=ЦШ,, 


Ф=ч/(И,-Ь,), 9=9.1,/И,, ОО (СЫ) 


где Х,У — безразмерные координаты, соответствующие координатам Х,У; [, — длина области 


И=ИиШ 


т! 


решения по оси; АИ — безразмерные скорости, соответствующие скоростям ИИ; И, — 


скорость потока на входе в воздушную камеру; Ф — безразмерный аналог функции тока; 
© — безразмерный аналог вектора вихря; С — безразмерная концентрация примеси; С. — 


концентрация источника; О — безразмерный аналог источника. 
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Для повышения сходимости решаемых уравнений в среде ЯехРОЕ-6.20 в модель был вве- 


дён поправочный коэффициент (/Ке ), который привёл исследуемые величины к одному поряд- 
ку. Система уравнений (1) приобрела вид: 


20) 2) : [20 29 
9х 09 Ве |0Х? ду? 


20 22 _ 
ь оыНИ С 











9Х>. 6’ 
с(0С) За 2 А 
1 С 6?С - 
—^+ Е . + |+^/Ае. 2 
9х 0У —5с.\Ве [5 = $ ы. 
И=д/ ах 
-У =дф/дУ, 


где И=0. /Ке, ИИ. /Ке, ф-Ф. Ве — скорости и функция тока системы уравнений с по- 


правочным коэффициентом (2). 
Граничные условия на входе в воздушную камеру: 


ФЕУ. /Ве, 9=0, С=0. 
Граничные условия на выходе: 
оф/ох =0, 9=0, 0б/оХ-У/..С.5с.ве. 


Граничные условия на стенках воздушной камеры [4, 5]: 


нижняя: Ф=0, и: оС/оХ =0; 

правая: ш-У./е, 9-=—^ __^, об/Х-=0; 

р = / 

верхняя: ш-У ./Ве, — _^ 0С/9Х=0; 
: УВе . Е? / 


левая: Ф=У,. /Ке, ЕР, 9С/5Х =0, 
\Ке .В 
где © — значение вектора завихрённости на границе; Ё — безразмерная длина отрезка от гра- 
ничной точки «0» до ближайшей к границе точки «1» (безразмерная величина приграничного 
слоя); Ч и Ч! — значения функции тока на граничной точке «0» и приграничной точке «1» 


соответственно. 
Численные решения на основе метода конечных элементов. В рассматриваемой модели 


присутствует параметр, связанный с формулировкой граничных условий (Ё). Были произведены 
расчёты определения диапазона изменения безразмерной величины приграничного слоя 


(В=В /Ё.), влияющего на сходимость решения. 
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На рис. 2 показаны результаты расчётов исследуемых параметров (функция тока, вектор 
завихрённости, скорости и концентрация) при значении приграничного слоя й =0,01.[, где 


[ =1 — безразмерная длина камеры. 
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Рис. 2. Результаты расчётов функции тока (а), вектора завихрённости (6), скорости (в) и концентрации (г) при значении 
приграничного слоя В =0,01-[ ‚где [=1 — безразмерная длина камеры 


На рис. 3 показаны результаты расчётов исследуемых параметров при значении погра- 
ничного слоя В =0,001.[. 

Таким образом, определена величина А, при которой система устойчива (рис. 2), и, как 
видно из рис. 3, а, уменьшение величины Й приводит к невозможности получить сходимость 
решения. 

В математической модели также менялась величина входной скорости. Диапазон значений 
скорости Ш» (от 0,1 м/с до 0,6 м/с) был принят в соответствии с санитарно-гигиеническими нор- 
мативами. Однако сходимость решения при данных параметрах скорости не наблюдалась. 

Сходимость наблюдалась в диапазоне скоростей (рис. 4), значительно ниже реальных ве- 
личин (Ш, =10“ м/си менее). 
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Рис. 3. Результаты расчётов во ЯехРОЕ функции тока (а), вектора завихрённости (6), скорости (в) и концентрации (г) при 


значении пограничного слоя й = 0,001.[ 
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Рис. 4. Результаты расчётов во НехРОЕ функции тока (а), вектора завихрённости (6), скорости (в) и концентрации (г) при 
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Заключение. 

1. Проведённые численные эксперименты доказали, что решение исследуемых уравнений 
в диапазоне нормированных скоростей неустойчиво. 

2. Необходим переход к более производительной программной среде, например АМ$У$ 
или ЗОПа МоГК$. 
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